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Abstract 

In these lectures, superconductivity in impure thin films close to the ab- 
solute zero of temperature is discussed. The behavior as function of the 
applied magnetic field and the amount of impurities suggests the presence 
of a superconductor-insulator transition at zero temperature. The theory of 
superconductivity in the limit where all the electrons become tightly bound 
in pairs is used to explain the main characteristics of the transition. In that 
limit, where the theory becomes equivalent to a phase-only theory, electron 
pairs exist on either side of the transition. The presentation is pedagogical 
in nature and includes exercises as a learning aid for those new to the field. 



PACS: 74.40.+k, 71.30.+h, 64.60.Fr 



1 Introduction 

The topic of these lectures is superconductivity in impure thin films close to the 
absolute zero of temperature. Such superconducting films may by increasing ei- 
ther the amount of impurities, or the applied magnetic field go over to an insulat- 
ing state. The transition is believed to signal the presence of a quantum critical 
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point at zero temperature [JI], [2p. Quantum critical phenomena differ from conven- 
tional critical phenomena taking place at finite temperature, in that quantum rather 
than thermal fluctuations are important in the critical region. Since quantum phase 
transitions occur at zero temperature, they cannot be induced by changing the tem- 
perature as can thermal phase transitions, and another parameter is to be varied to 
trigger the transition. It also means that one can never tune right through the tran- 
sition as experiments are necessarily done at finite temperature. For low enough 
temperature though, the presence of a quantum critical point can nevertheless be 
detected by using finite-size scaling. 

A superconducting and an insulating state seems an unlikely combination to 
be present in the same system. Whereas superconductivity needs an attractive 
interaction for the pairing between electrons, the insulating state needs a repulsive 
interaction. It is a priori not clear how these two requirements can be fulfilled in a 
single system. Adding to the bewilderment is the striking similarity in the current- 
voltage characteristics in both phases close to the transition. By interchanging the 
current and voltage axes in one phase, a characteristic obtained at some value of 
the applied magnetic field say, can be mapped onto a characteristic of the other 
phase obtained at a different value of the field. This means that, although the 
physical properties of the superconducting and the insulating state are completely 
different, there is a close connection between the conduction mechanisms in the 
two phases. A last defining aspect of the superconductor-insulator transition is 
the presence of a 1/r Coulomb potential, which at low charge-carrier densities 
becomes very strong. 

In these lectures, we will argue that the main aspects of the superconductor- 
insulator transition just mentioned can be accounted for by a single theory. Namely, 
the theory of superconductivity in the limit where all the electrons become tightly 
bound in pairs — the so-called composite boson limit [|3], |J]. This is the opposite 
limit of the conventional weak-coupling BCS limit, where only electrons in a thin 
shell around the Fermi surface become loosely bound in Cooper pairs. In the com- 
posite boson limit, the superconducting state displays a quantum phase transition 
to an insulating state characterized not by an unbinding of electron pairs, but rather 
by a quenching of the condensate of composite bosons. In other words, electron 
pairs exist on both sides of the transition. In the superconducting state they are 
Bose-Einstein condensed, while in the insulating state they are localized. The at- 
tractive interaction responsible for the binding of the electrons in pairs translates 
into a repulsive interaction between the composite bosons. The theory describ- 
ing these bosons is the Bogoliubov theory of superfluidity, which is equivalent 
to a phase-only theory. This implies that the superconductor-insulator transition 
can be described by a phase-only effective theory without ignoring any degrees of 
freedom [§, |p . 

To be fair, we should mention at this point that no consensus exists to what 
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extend the phase-only theory can be applied to the superconductor-insulator tran- 
sition. Further experimental and theoretical studies are required to settle this point. 

1.1 General Scaling 

Since quantum rather than thermal fluctuations are relevant in a quantum phase 
transition, one has to work in spacetime rather than just in space as is appropriate 
for thermal phase transitions in equilibrium. As a result, in addition to a diverging 
correlation length £, quantum phase transitions have also a diverging correlation 
time £ t . They indicate, respectively, the distance and time period over which the 
system fluctuates coherently. The two are related, with the diverging correlation 
time scaling with the diverging correlation length as 

6 ~ f , (i) 

where z is the so-called dynamic exponent. It is a measure for the asymmetry 
between the time and space directions close to the critical point. The dynamic ex- 
ponent is to be added to the set of critical exponents used to characterize a thermal 
phase transition. Since such transitions have only two independent exponents, a 
quantum phase transition is specified by three independent exponents. 

The traditional scaling theory of thermal, continuous phase transitions in equi- 
librium, first put forward by Widom [J7p, is easily extended to include the time di- 
mension [|8[] because the relation ([]]) implies the presence of only one independent 
diverging scale. 

Let S = K—K c , with K the parameter triggering the phase transition, measure 
the distance from the critical value K c . A physical observable O at the absolute 
zero of temperature depends on K as well as on other variable, such as an external 
field, energy, or momentum. Let us denote these other variables collectively by V. 
In the critical region close to the critical point, O can be written as 

0(T,K) = ^0(t), (T = 0), (2) 

where do is the scaling dimension of the observable O, £ ~ \5\~ u , with v the 
correlation length exponent, and T is obtained from T by rescaling it with factors 
of the correlation length, so that T is independent of that scale. To be specific, if an 
external field scales as T ~ £ dr , then the rescaled field is defined as T = £~ dr T. 
The right side of Eq. (g) does not depend explicitly on K, but only implicitly 
through £. 

The data of the observable O as function of an external field T obtained at 
different values of the parameter K triggering the transition can be collapsed on 
a single curve when instead of 0(T, K), the rescaled quantity \8\ ud °0(T, K) is 
plotted as function not of T, but of T. Indeed, because of Eq. (g), the combination 
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I^I^C^r, K) depends only on f and is thus independent of the distance from 
the critical point. This is used to determine critical exponents experimentally. By 
rescaling the vertical and horizontal axis of the plot with \6\ x and \5\ y ', respectively, 
the best collapse obtained at some value x and y give the combination of critical 
exponents vd Q = yo and ud r = x . 

Since a physical system is always at some finite temperature, we have to in- 
vestigate how the scaling law @ changes when the temperature becomes nonzero. 
The easiest way to include the temperature in a quantum theory is to go over to 
imaginary time r = it, with r restricted to the interval < r < 1/T. The tempo- 
ral dimension thus becomes of finite extend. The behavior at a finite temperature 
is still controlled by the quantum critical point, provided the correlation time sat- 
isfies < 1/T. If this condition is fulfilled, the system will not see the finite 
extent of the time dimension. This is what makes quantum phase transitions ex- 
perimentally accessible. Instead of the zero-temperature scaling (0), we now have 
the finite- size scaling 

0(T, K, T) = T- d °/ z O(T T , &T), (TV 0), (3) 

where instead of using the correlation length to convert dimensionfull quantities 
in dimensionless ones, the temperature is used: fY = T dr ^ z T. 



Notation 

We adopt Feynman's notation and denote a spacetime point by x = x M = (t, x), 
fi = 0, 1, • • • , d, with d the number of space dimensions, while the energy k 
and momentum k of a particle will be denoted by k = = (k , k). The time 
derivative do = d/dt and the gradient V are sometimes combined in a single 
vector <9 M = (do, — V). The tilde on <9 M is to alert the reader for the minus sign 
appearing in the spatial components of this vector. We define the scalar product 
k ■ x = k^x^ = k t — k • x and use Einstein's summation convention. Because of 
the minus sign in the definition of the vector it follows that d^a^ = <9 a + V • a, 
with an arbitrary vector. 

Integrals over spacetime are denoted by 




dtd d x, 



while those over energy and momentum by 

f _ f f d ^o d d k 
Ik ~ A ,k ~ J 2tt (2vr) d ' 
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When no integration limits are indicated, the integrals are assumed to run over all 
possible values of the integration variables. Similarly, for functional integrals we 
use the notation 

We will work in natural units with the speed of light, Boltzmann's constant 
&b, and Planck's constant % set to unity. 

These lectures include exercises, which are clearly marked. Most of the solu- 
tions can be found in Ref. [Q]. 



2 Superconductivity 

In this section we study the theory of superconductivity in the limit where all the 
electrons become tightly bound in pairs [|], f|]. The composite boson limit is to 
be distinguished from the usual weak-coupling BCS limit, where only electrons 
(of opposite momentum) in a thin shell around the Fermi surface become loosely 
bound in Cooper pairs. 



2.1 BCS Theory 

As starting point we take the microscopic BCS model specified by the Lagrangian 
PD 



c = v T *H-eHv)]v T +^[^o-e(-iv)]^-A^^^V'T 

:= C + C- mt , (4) 

where C is the free theory, and C mt = — A^? i^l ip^ ipi is a local electron-electron 
interaction term, representing the effective, phonon mediated, attraction between 
electrons with coupling constant A < 0. The field is an anticommuting 
field describing the electrons with mass m and spin up (down), while £(— zV) = 
e(— zV) — fJL, with e(— z'V) = — V 2 /2m, is the kinetic energy operator with the 
chemical potential [i characterizing the ensemble of fermions subtracted. The 
theory is invariant under global U(l) transformations under which 

^ -> e"V CT , (5) 



with a =1,1 and a sl constant transformation parameter. The superconducting 
state is characterized by a spontaneous breakdown of this symmetry. 

To investigate the superconducting state, we integrate out the fermionic de- 
grees of freedom at the expense of a newly introduced auxiliary field, which is 
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better equipped to describe this condensed matter state. It is convenient to first 
introduce the two-component field 



In this so-called Nambu representation, £ becomes 

where the anticommuting character of the electron fields is used and total deriva- 
tives are neglected. We next replace the electron-electron interaction with an ex- 
pression involving the auxiliary field A 

VI V»i V>T - A* Vt + V T * V! A - j | A| 2 . (8) 

The original interaction is regained when the auxiliary field is integrated out. 
Physically, A, representing a product of two electron fields, describes electron 
pairs. We shall therefore refer to it as pair field. With this replacement, the parti- 
tion function 

Z = ( exp (z / c), (9) 

becomes 

Z = / exp ( | / |A| 2 ) (10) 



x exp 



, t I ido - e(-<V) -A \ . 
V I -A* zdo + £(-iV i ^ 



Since the fermion fields appear quadratically now, they may be integrated out to 
yield an effective action S e s for A and A* 

^[a-.ah -mm po 4 (p) ), (ID 

where p = id and £(p) = e(p) — ji, with e(p) = p 2 /2m and p = — iV, so that 
the partition function becomes 

Z = |^ A exp^ cff [A*,A] + ^|jA| 2 ). (12) 



The trace Tr appearing in Eq. < Jll[ > denotes the trace over discrete indices as well 
as the integration over spacetime and over energy k and momentum k. 
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In the mean-field approximation, the functional integral over the pair field 
in Eq. ( |12| ) is approximated by the saddle point, where only the extremal value, 
satisfying the equation 

~F~\~*7 \ = -T A 0[X), (13) 

oAq(x) A 

is included. For a system homogeneous in spacetime, the pair field is a constant 
A, and Eq. (|T3l) reduces, after passing to the Fourier representation, to the BCS 
gap equation JTH : 

l r l 



A " Jkk%-E 2 (k)+iri 

(14) 



2 A E(k) ' 

Here, i] is an infinitesimal positive constant which is to be set to zero at the end of 
the calculation, and 



E(k) = ^ 2 (k) + |A | 2 (15) 

is the spectrumQ of the fermionic excitations. A nontrivial solution to the gap 
equation signals the spontaneous symmetry breakdown of the global U(l) sym- 
metry (|5]). 

2.2 Composite Boson Limit 

For studying the composite boson limit, it proves prudent to swap the coupling 
constant A in the gap equation dl4| ) for a more convenient parameter, namely the 
binding energy e a of an electron pair in vacuum [JTTp. Both parameters character- 
ize the strength of the electron-electron interaction. To establish the connection 
between the two, let us consider the Schrodinger equation for the problem at hand. 
In reduced coordinates, it reads 



V 2 

— + A<5(x) 
m 



^(x) = -e a , (16) 



where the reduced mass is m/2 and the delta-function potential, with A < 0, 
represents the attractive local electron-electron interaction £ int in Eq. (Q). We 
stress that this is a two-particle problem in vacuum and not the famous Cooper 
problem of two interacting fermions on top of a filled Fermi sea. The equation is 
most easily solved in the Fourier representation, yielding 

= -kv^' (0) - (17) 



'To avoid confusion, let us repeat that the bar in Ao indicates that the pair field is a constant, 
while the subscript indicates that it satisfies the extremal condition dl3b. 
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or 

1 r 1 

-T = / V27 I • (18) 

A Jk k z /m + e a 

This bound-state equation allows us to replace the coupling constant A with the 
binding energy e a . When substituted in the gap equation (|14|), the latter becomes 

Jkk 2 /m + e a ~ 2Jk£(k)' (1 } 

By inspection, it follows that this equation has a solution given by [|3], |]] 

A -> 0, // -> -|e B , (20) 

where it should be noted that, in contrast to the weak-coupling limit, the chemical 
potential characterizing the ensemble of fermions is negative here. This is the 
composite boson limit. 

To appreciate the physical significance of the specific value found for the 
chemical potential in this limit and also its name, observe that the spectrum E B (q) 
of the two-fermion bound state measured relative to the pair chemical potential 2/x 
reads 2 

£ B (q) = -r--/i B , (2D 

Am 

where lib is defined as /i B = e a + 2/x and may be understood as the chemical 
potential characterizing the ensemble of composite bosons. The negative value 
for jx found in Eq. ( |20| ) is precisely the condition for a Bose-Einstein condensation 
of an ideal gas of composite bosons in the q = state. 

Including quadratic terms in A , we obtain as solution to Eq. ( |19| ) 

^=~ea + (l-d/4)i^L. (22) 
This leads to the chemical potential 

IA I 2 

^ B =(2-d/2)i-^-, (23) 
characterizing the now interacting Bose gas. 



2.3 Renormalization 

For a system homogeneous in spacetime, so that the field Ao(x) is constant, the 
effective action < JTT| ) is readily evaluated. Disassembling the argument of the log- 
arithm as 



Po - UP) ~ A o \ Po 






(24) 
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and expanding the second logarithm in a Taylor series, we recognize the form 

^a ; ,a ,= ~^M ip) P0+ ° f(p) ) 



I A, 
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-zTrln 1- 2 ' , (25) 

apart from an irrelevant constant. Again passing to the Fourier representation, 
and carrying out the integral over the loop energy k , we obtain the effective La- 
grangian 

C cS = jf [E(k) - f (k)] . (26) 

Exercise: Derive this result from Eq. (p5|), using contour integration. Rather than integrating the 
logarithms in that equation, one better first differentiate it with respect to the chemical potential 
as the integral over the loop energy fco becomes easier that way. In the end one integrates the 
resulting expression again over the chemical potential to obtain the desired result I 



To the one-loop result (|2~6]), the tree term |A | 2 /A is to be added. Expanding 
E(k) in Eq. ( |26l ) in a Taylor series, we see that the effective Lagrangian also 
contains a term quadratic in A . This term amounts to a renormalization of the 
coupling constant. Specifically, the renormalized coupling constant A r is to this 
order given by 

1 1 1 

v + ^ 



A r A 2 A f (k) 

1 T(l - d/2) m d / 2 

A + (4tt)<V 2 e i" d / 2 ' V n 

with, as is appropriate in the composite boson limit, 

£(k) = e (k) + |e a , (28) 

and dimensional regularization is used in evaluating the momentum integral. Be- 
cause of Eq. dT8|), which may be viewed as the defining equation of the parameter 
e a , the right side of Eq. ( |2"7| ) is zero. This implies that in the composite boson 
limit, A r — > — oo so that we indeed have tightly bound pairs here. 



3 Bogoliubov theory 

In this section we show that in the composite boson limit, the effective theory 
obtained after integrating out the fermionic degrees of freedom, is the Bogoliubov 
theory of superfluidity in an interacting Bose gas. The gas consists of composite 
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bosons with a mass twice the electron mass. The system is known to undergo 
a quantum phase transition from the superfiuid to a (Mott) insulating state JT2|]. 
When translated back to the fermionic theory, this transition corresponds to one 
where the condensate is drained of composite bosons, without breaking them up. 
We include impurities in the Bogoliubov theory to show that it leads to localization 
without destroying the superfiuid state completely. The insulating state, which is 
now a result not only of repulsive interactions, as in a Mott insulator, but also of 
(Anderson) localization, is called an Anderson-Mott insulator. 



3.1 Derivative Expansion 

We next wish to relax the assumption of homogeneity in spacetime and consider 
a spacetime-dependent pair field. To this end, we study the effective action ( [TT| ) 
and expand A(x) around the constant value A satisfying the gap equation (|R|), 



A(x) = A + A{x). 



(29) 



We obtain in this way, 

See 

where Go is the correlation function, 

G (k) = 



OO 



G (p) 



A 

A* 



(30) 



k - £(k) -A \ 

-a* h+av J 

1 ( k e ikor i_ + ^(k) 

A* 



(31) 



A 

k e~ ik ^ - 



kl - £ 2 (k) +i v \ A* k e^" - £(k) J ' 

The exponential functions in the diagonal elements of the correlation function are 
additional convergence factors needed in nonrelativistic theories [|13p. 

When evaluating the effective action (^0|), the precise meaning of the trace Tr 
appearing there should be kept in mind. Explicitly, it is defined as 



S e f{ = —iTr In [K(p, x)] = — itrln K(p,x)5(x — y)\ 



\y=x 



(32) 



where the trace tr is the usual one over discrete indices. We abbreviated the matrix 
appearing in Eq. < JTT] ) by K(p, x) so as to cover the entire class of actions of the 
form 



S 



ip'(x)K(p, x)ip(x) 



(33) 
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The delta function in Eq. arises because K (p, x) is obtained as a second 
functional derivative of the action 

WWW) = k <p> x W x -v)\v=.> ^ 

each of which gives a delta function. Since the action contains only one inte- 
gral over spacetime, one delta function remains. Because it is diagonal, the delta 
function may be taken out of the logarithm and Eq. d32[ ) can be written as 

y=x 

= -itr / I e lk ' x \n[K(p,x)}e- ik - x . (35) 

Jx Jk 

In the last step, we used the integral representation of the delta function: 

S(x) = I e~ ik - x , (36) 



^ = -itr / In [K(p, x)} 8{x-y)\ 

Jx 



shifted the exponential function exp(ik ■ y) to the left, which is justified because 
the derivative p^ does not operate on it, and, finally, set equal to x^. We thus 
see that the trace Tr in Eq. (TJ2|) stands for the trace over discrete indices as well 
as the integration over spacetime and over energy and momentum. The integral J k 
arises because the effective action calculated here is a one-loop result with the 
loop energy and momentum. 

The integrals in Eq. (|35|) cannot in general be evaluated in closed form because 
the logarithm contains energy-momentum operators and spacetime-dependent func- 
tions in a mixed order. To disentangle the integrals resort has to be taken to a 



derivative expansion [ [T4| ] in which the logarithm is expanded in a Taylor series. 
Each term contains powers of the energy-momentum operator p^ which acts on 
every spacetime-dependent function to its right. All these operators are shifted to 
the left by repeatedly applying the identity 

f(x)p ll g(x) = (jv - iB ll )f(x)g(x), (37) 

where f(x) and g(x) are arbitrary functions of spacetime and the derivative <9 M = 
(do, — V) acts only on the next object to the right. One then integrates by parts, so 
that all the p M 's act to the left where only a factor exp(ik ■ x) stands. Ignoring total 
derivatives and taking into account the minus signs that arise when integrating by 
parts, one sees that all occurrences of p^ (an operator) are replaced with k^ (an 
integration variable). The exponential function exp(— ik ■ x) can at this stage be 
moved to the left where it is annihilated by the function exp(ik ■ x). The energy- 
momentum integration can now in principle be carried out and the effective action 
be cast in the form of an integral over a local density £ eff : 

See = / C c s- (38) 
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This is in a nutshell how the derivative expansion works. 



Exercise: Apply the derivative expansion to the Lagrangian 

£=\{d^) 2 -|mV -\\<P\ (39) 

with <fi a real scalar field. The theory has a Z2 symmetry under which the scalar field flips sign: 
4> — > 4>' — — <fi. Show that the effective theory in two space dimensions is given by 

C cS = \Z($)(d^f - V ff(^), (40) 



with 



VeffW = --^(m 2 + iA^) 3 / 2 , 

iZTT 

1 \ 2 6 2 

= 1927T ( TO 2 + l A02) 3/2- < 41 > 



3.2 Map onto Bogoliubov Theory 

We are interested in terms in the effective action ( [301 ) quadratic in the field A. 
Using the derivative expansion, we find 

6ff [Q) 2 pi - ^(p) (po + g ) 2 " £ 2 (P " q) ^ } 



{a 2 a* a* + [p + e(p)][pd + go - e(p - q)] aa* 

+A* 2 AA + [po - £(p)][p + go + f (p - q)]A*A}, 



where g M = id^. Let us first ignore the derivatives in this expression. After 
carrying out the integral over the loop energy k and using the gap equation ([R]), 
we obtain 

£ (2) (0) = -~ I ^) (A 2 A* 2 + Af A 2 + 2| A | 2 | A| 2 ) . (43) 

In the composite boson limit, where the spectrum of the fermionic excitations 
is given by Eq. ([28]), the integral over the loop momentum becomes elementary, 
yielding 

1 _4r(3-d/2) rod/2e , /2 _ 3 _ 



E 3 {k) (47T) d / 2 

We next consider the terms in Eq. ( |42~| ) involving derivatives. Following Ref. 
[ [131 ] we set A to zero here. The integral over the loop energy is easily carried 
out, with the result 

C (2) (q) = -- f - AA* 

2 A q - k 2 /m + 2/i - q 2 /4m 

If 1 ~ ~ 

— / A*A. (45) 

2 A —go — k 2 /m + 2/i — q 2 /4m 
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In the composite boson limit, the remaining momentum integrals become elemen- 
tary again and after expanding in derivatives we find 



L 



k q — k 2 /m — e a — q 2 /4m 



' (46) 



r(l - d/2), /2 , /2 _ x r(2 - d/2), nM / q 2 



_ md ,2 ed ,2-l _ ZSZ J/^ /2^/2-2 % 



( 4vr )d/2 - ( 47r )d/2 a yiu ^ j ■ 

When combined with the tree term |A| 2 /A, the first term at the right side of this 
equation yields the renormalization ( |2"7| ) of the coupling constant. The second 



term at the right side of Eq. ( |46j ) gives, when combined with the contribution ([43]), 
the result [fBt], 



£ (2, = i £ |^". i/ M /2 - 2 * , M„( 9 )*, *=(f.), (47) 

where M (q) stands for the 2x2 matrix, 

M (q) = 

f q - q 2 /4m - (2 - d/2)\A Q \ 2 /e a -(2 - d/2)A 2 /e a \ 

{ -(2 - d/2)Af/e a -q - q 2 /4m - (2 - d/2)|A | 2 Aa J ' 

(48) 

As we shall show now, this is the Bogoliubov theory of superfluidity in an in- 
teracting Bose gas. That is to say, after integrating out the fermionic degrees of 
freedom from the theory of superconductivity, we obtain in the composite boson 
limit a theory describing a gas of repulsively interacting (composite) bosons. 



3.3 Quantum Phase Transition 



The Bogoliubov theory is specified by the Lagrangian [ |T6| ] 



C = (j)*[id - e(-iV) + ^ B ]0 - A B |0| 4 , (49) 

where /xb is the chemical potential characterizing the Bose gas. The self-coupling 
is taken to be positive, Ab > 0, so that the local interaction is repulsive. 

At the mean-field, or classical level, where quantum fluctuations are ignored, 
the theory (^) undergoes a phase transition when the chemical potential changes 
sign. For /j B > 0, the global U(l) symmetry is spontaneously broken by a nontriv- 
ial ground state, while for /i B < the symmetry is unbroken. The change in /i B 
can be induced by varying the temperature, as in a thermal phase transition, but it 
can also be induced at zero temperature by varying, for example, the number of 
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charge carriers, or the amount of impurities. The zero-temperature quantum phase 
transition describes a transition between a superfiuid and an insulating state [fi~2|]. 

The ground state of a system homogeneous in spacetime is obtained by con- 
sidering the shape of the potential energy 

V = -/i B |0| 2 + A B |0| 4 . (50) 
For /i B > it indeed has a minimum away from the origin = given by 

I J 12 1 A*B 



^0 



2A 



(51) 



B 



and the potential becomes 



(52) 

Since the total particle number density n B is represented by 

n B (x) = \<f ) (x)\ 2 , (53) 

the quantity n : = |0 O | 2 physically denotes the number density of particles re- 
siding in the ground state. A nonzero value for n thus signals Bose-Einstein 
condensation. For a homogeneous system in its ground state, we see that at the 
mean-field level no = n so that all the particles reside in the condensate. This will 
change when quantum fluctuations are included as a result of which particles are 
knocked out of the condensate (see below). 

To account for the nontrivial ground state, we introduce the shifted fieldQ (j)(x) : 

(j>{x) = Q + 0(x). (54) 



The terms in the Lagrangian (|49|) quadratic in this shifted field may be cast in the 
matrix form 

C = \¥M Q {p% <§ = d \ (55) 



with 



M (p) 



Po - e(p) + u B -4A B |0o| — 2A B 0o 



v -2\ B (pf -p o -e(p)+/i B -4A B |0 o | 2 J ' 

(56) 

Taking into account only the quadratic terms in the field and neglecting higher- 
order terms, as we just did, is known as the Bogoliubov approximation. 

2 Similar as before, the bar in <fi denotes a constant value of the field, while the subscript 
indicates that it satisfies the mean-field equation (pl|). 
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Comparing this expression with Eq. ( fi8| ) obtained in the composite boson limit 
after integrating out the fermionic degrees of freedom from the theory of super- 
conductivity, we conclude that the composite bosons have — as expected — a mass 
m-Q = 2m twice the fermion mass m, and a small chemical potential given by 



Eq. (p3|), which we there derived from the gap equation. It also follows that the 



number density of composite bosons condensed in the ground state reads 
and that the interaction Ab between the composite bosons is 
or, using Eq. (JT§), 

Note that the parameter A(< 0) characterizing the attractive electron-electron in- 
teraction appears below d = 2 with a minus sign here, leading to a repulsive 
interaction between the composite bosons. (In the next subsection, we will see 
that d = 2 is the upper critical dimension of the T = Bogoliubov theory.) This 
brings us to the important conclusion that for d < 2 the same interaction respon- 
sible for the formation of electron pairs, is also responsible for the stability of the 
superfluid state, and when this state ceases to exist, for that of the insulating state, 
which both need a repulsive interaction. 

The quantum phase transition encoded in the Bogoliubov theory corresponds, 
when translated back to the fermionic theory, to one where the condensate is 
drained of composite bosons, without breaking them up. In other words, com- 
posite bosons exist on both sides of the transition, either condensed (superfluid 
state) or localized (insulating state) [EL Bp. 



3.4 Beyond mean-field theory 

We can continue now and improve on the usual mean-field approximation of the 
theory of superconductivity, where the functional integral over the pair field in the 
partition function ( TT2| ) is approximated by the saddle point, by integrating out the 
field ^ in Eq. (p7|), or to simplify notation, the field $ in Eq. (|55|). This leads to 
the effective potential 

V eff = -^tr / ln[M (A;)] = \ [ E(k). (60) 
2 Jk 2 A 
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Here, E(k) is the famous single-particle Bogoliubov spectrum [|T5p, 



E(k) = v/e 2 (k)+2/i B e(k) 



/ e 2 (k)+4A B |0o| 2 e(k). (61) 
In the limit of large momentum, the spectrum behaves in a way 

E(k) ~e(k) + 2A B |0o| 2 (62) 

typical for a nonrelativistic particle of mass m moving in a background medium, 
provided by the condensate in this case. The most notable feature of the Bogoli- 
ubov spectrum is that it is gapless, behaving for small momentum as 

£7(k)~c|k|, (63) 

with c = \J pl-q/tu. 

Exercise: Carry out the integral over the loop energy ko in Eq. ( |60| ) using contour integration and 
show that this leads to the right side of that equation. This is best done by first differentiating the 
expression with respect to the chemical potential /jb, and in the end integrating the result again 
with respect to /xb . 



The integral over the loop momentum in Eq. ( |60| ) can be carried out using the 
integral representation of the Gamma function 

1 1 r°° dr 

T z e~ aT . (64) 



a z T(z) Jo t 

In arbitrary space dimension d this yields, using dimensional regularization: 

V - T m*/*u d '™ T - r(l-<*/2)lW2 + l/2) _ 
V cS --L d m ^ , U- 2nd /2 + i/2 Tid/2 + 2) ■ ( 65 > 

For d — 2, the effective potential diverges. To investigate this, we set d = 2 — e, 
with e small and positive, and expand V e g around d — 2, giving 

Vcff = "4^^' (66) 

with k an arbitrary renormalization group scale parameter which enters for di- 
mensional reasons. If the Bogoliubov spectrum had not been gapless, but had an 
energy gap instead, this parameter would have appeared in Eq. ((6^) in the place 
of k. As always in dimensional regularization, the divergence shows up as a pole 
in e. Comparing the one-loop contribution with the classical contribution (|52|), 
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we conclude that Eq. (153) leads to a renormalization of the coupling constant A, 



yielding the renormalized coupling A r QT7| ] 



1 1 ml 

T = j + , (67) 

A r A 7T e 

where A = \/k £ / 2 and a similar definition for A r . The quantum critical point is 
approached by letting the renormalized group scale parameter k — ► 0. For fixed 
coupling A, it then follows that upon approaching the critical point, the renormal- 
ized coupling tends to A r — > ixe/m. For d < 2, or equivalently e > 0, the fixed 
point is nontrivial. In the limit d — > 2, A r — > and the theory becomes Gaussian, 
identifying d = 2 as the upper critical dimension. 

Due to quantum fluctuations not all the particles are known to reside in the 
condensate [fTgp. Specifically, in d space dimensions, the (constant) particle num- 
ber density n at the one-loop level is given by [ |19| ] 



n = |0 O | 2 - 2 rf/2 - 2 ^W /2 Af |0o| d (68) 

Since the quantum-induced term is positive for 1 < d < 4, the number of particles 
residing in the condensate given n is reduced compared to the classical result 
n = |0o| 2 - This shows that due to quantum fluctuations, particles are knocked out 
of the condensate. 

Exercise: Derive Eq. (^8j). In doing so, one should not use the mean-field equation (|T]) too early, 
and instead work with the more general single-particle spectrum 



E(k) = yj [e(k) - fi Q + 4A O |0| 2 ] - 4A^|4 . (69) 

It reduces to the Bogoliubov spectrum when the mean-field equation is used. 

Despite that due to quantum fluctuations not all the particles reside in the con- 
densate, all the particles do in the absence of impurities and at zero temperature 
participate in the superflow, and move on the average with the superfluid velocity. 
Put differently, the superfluid mass density p s is given by the total particle number 
density n: p s = mn. 

Exercise: Prove this. To this end, assume that the entire system moves with a velocity u relative to 
the laboratory system. As in standard hydrodynamics, the time derivative in the frame following 
the motion of the fluid is 8q + u • V. Also assume that the condensate moves with the superfluid 
velocity v s and boost the field: 

tf>{x) -» 4t(x) = e' imVs - x <Hx). (70) 

Show that when incorporated in the Lagrangian (^) of the interacting Bose gas, these two changes 
result in a change of the chemical potential 

Mb -> Mcff := Mb - \mv a ■ (v s - 2u). (71) 
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Show that the resulting Bogoliubov spectrum and thermodynamic potential are given by the pre- 
vious results ( |6l| ) and ( |65| ) with this replacement. 

The momentum density, or equivalently, the mass current g of the system is obtained in this 
approximation by differentiating the effective potential with respect to u. Show that 

g = p s v s , (72) 

with p s — mn the superfiuid mass density. 



3.5 Impurities 

One of the ways to trigger a superconductor-insulator transition is to change the 
amount of impurities. This means that, e.g., the correlation length £ diverges as 
| a* — a\~ u when the parameter a characterizing the impurities approaches the 
critical value a*. 

To account for impurities, we include a term 

£ a = ^(x)|0(:r)| 2 (73) 

in the bosonic theory (ff9[), where if)(x) is a space-dependent random field with a 
Gaussian distribution 



P[ip] = exp 



-- / ^ 



a Jx 



(74) 



characterized by the impurity strength a (a alluded to above is a rescaled version 
of a). Notice that the random field does not depend on time. This is because it 
is introduced to mimic impurities, which are randomly distributed in space, not in 
time. 

We shall treat the impurities in the so-called quenched approximation [|8|], 
where the average of an observable 0(0*, 0) is obtained as follows 

(0(0*, 0))= f PM(O(0*,0)>^ (75) 

with (0(0*, 4>))^ indicating the grand-canonical average for a given impurity con- 
figuration. That is to say, first the ensemble average is taken for fixed -0, and only 
after that the averaging over the random field is carried out. 

In terms of the shifted field (Q), the random term ( f73"| ) becomes 

C a = V(x)(|0 o | 2 + |0| 2 + 000* + 000). (76) 

The first two terms lead to an irrelevant change in the chemical potential, so that 
only the last two terms need to be considered, which can be cast in the matrix 
form 

C a = V^(x) $ $, $ = f t) ■ (77) 
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The integral over <3> is Gaussian in the Bogoliubov approximation and there- 
fore easily performed to yield an additional term to the effective action 

Sa = ~[ H*)HGo(x-y)%Hy), (78) 

Z J x,y 

where the correlation function Go is the inverse of the matrix M introduced in 
Eq. (|56|). To proceed, we pass to the Fourier representation: 



Go(x-y) = f e- lk -^G (k) (79) 
^(x) = ^ e iko V(k). (80) 
The contribution to the effective action then appears in the form 

Sa = ~~ I |^(k)| 2 <5 G(0,k)<i> . (81) 
2 A 

Since the random field is Gaussian distributed, the average over this field repre- 
senting quenched impurities yields: 

(|^(k)| 2 ) = |fi«, (82) 

with n the volume of the system. The remaining integral over the loop momen- 
tum in Eq. < ]8T] ) is readily carried out to yield in arbitrary space dimensions the 
contribution to the Lagrangian 

(C a ) = -r(l - d/2) (-) |0o| 2 (6A B |0o| 2 - ^-'a. (83) 

The divergence in the limit d — > 2 shows that also in the presence of impurities, 
the two-dimensional case is special. This expression can be used to obtain the 
additional depletion due to impurities. To this end, we differentiate it with respect 
to the chemical potential, giving 



_ d(C a ) 2 rf / 2 - 5 r(2 - d/2) d/2 . d /2-2_ d/2 -i faA . 
n a = — = jj- m 1 A B ' n ' a, (84) 

where we recall that n = |0 O | 2 denotes the (constant) number density of particles 
residing in the condensate. Because this contribution is positive, it amounts to an 
additional depletion of the condensate. The divergence in the limit Ab — > for 
d < 4 signals the collapse of the system with impurities when the interparticle 
repulsion is removed. 

To determine the superfluid mass density p s in the presence of impurities, we 
replace, as in the last exercise of Sec. |3T4"| , // B with fi cS defined in Eq. ( J7T| ) and 
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ido with ido — (u — v s ) ■ (— iV) in the contribution ( |8~T| ) to the effective action. 
Differentiating it with respect to the externally imposed velocity, — u, we find to 
linear order in the difference u — v s : 

g = p s v s + p n u, (85) 

with the superfluid and normal mass density pT) ] 

/_ 4_ \ _ 4 
p s = m [n — —n a ] , p n = ~mn a . (86) 



As expected, p s ^ ran in the presence of impurities. Moreover, the normal mass 
density is a factor A/d larger than the mass density mn a knocked out of the con- 
densate by the impurities. For d = 3 this gives the factor | first found in Ref. []2~2|]. 



As argued there, this indicates that part of the zero-momentum state belongs not 
to the condensate, but to the normal fluid. Being trapped by the impurities, the 
fraction (4 — d)/d x n a of the zero-momentum state is localized. 

This is an important conclusion as it shows that the phenomenon of localiza- 
tion can be accounted for in the Bogoliubov theory of superfluidity by including 
a random field, without necessarily destroying that state. 



4 Phase-only theory 

In this section we show that the Bogoliubov theory, which we obtained in the 
composite boson limit after integrating out the fermionic degrees of freedom from 
the theory of superconductivity, contains only one degree of freedom, viz. the 
phase of the order parameter. Physically, it describes the Goldstone mode of the 
spontaneously broken global U(l) symmetry. In the context of superconductivity, 
this mode is called Anderson-Bogoliubov mode. The Bogoliubov theory may 
therefore, at least in the superfluid state, be represented by a phase-only effective 
theory. We continue to account for the 1/r Coulomb potential in the effective 
theory and give general scaling arguments for the physical quantities represented 
by that theory. 



4.1 Derivation 



It was first shown by Beliaev Q23Q that the gaplessness of the single-particle spec- 
trum first found by Bogoliubov at the classical level persists at the one-loop order 
and later proven by Hugenholtz and Pines [ [24] ] to hold to all orders in perturba- 
tion theory. In fact, as was proven by Gavoret and Nozieres P5|], the Bogoliubov 
spectrum is identical to that of the Goldstone mode accompanying the sponta- 
neous breakdown of the global U(l) symmetry, thus explaining its gaplessness. 
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Also from the perspective of degrees of freedom, this conclusion makes sense. 
Although the normal phase is described by a complex 0-field, having two com- 
ponents, it contains only one degree of freedom [E5p. This is because the energy 
E(k) ~ k 2 is always positive. As a result, only positive energies appear in the 
Fourier decomposition of the field, and one needs — as is well known from stan- 
dard quantum mechanics — a complex field to describe a single spinless particle. 
In the superfiuid phase, on the other hand, where E 2 (k) ~ k 2 , the Fourier de- 
composition contains positive as well as negative energies so that a single real 
field suffices to describe this mode. In other words, although the number of fields 
is different, the number of degrees of freedom is the same in both phases. This 
implies that the superfiuid state can be described by a phase-only theory as it cap- 
tures all the degrees of freedom, ignoring vortices for the moment which are easily 
incorporated in the theory as will be discussed in the next section. 

To obtain the phase-only theory, we set, cf. Eq. (|53|) 

<f>(x)=e i *W(fa + $), (87) 

with tp(x) a background field representing the Goldstone mode accompanying the 
spontaneous symmetry breakdown of the global U(l) symmetry. Inserting this in 



the Lagrangian ( j49D and expanding it, we obtain 

£ (2) = -V - \fo\ 2 U - U(fa4>* + 0*0) - \ B \M 2 (M* + (88) 
where the field U (x) stands for the combination 

U(x) = d ip(x) + ^-[V V (x)] 2 . (89) 

In deriving Eq. (|88l), we used the mean-field equation /zb = 2Ab|0o| 2 - We con- 
tinue to integrate out the tilde field (which is tantamount to substituting its field 
equation back into the Lagrangian) to obtain the phase-only theory 

£ eff = -nU(x) + ]u(x)—U(x), (90) 

4 Ab 

where we ignored the irrelevant constant term Vo and substituted |0 O | 2 (= ^o) = n 
to this order. Using the mean-field equation again, we can write the coefficient of 
the last term as: 

1 1 1 n 1 9 

-— = 7 = -n 2 n, (91) 

4 A B 2 mc 2 2 v 

with c the speed of sound introduced in Eq. (Q). Standard thermodynamics relates 
c to the compressibility k via 

K = — (92) 
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The phase-only theory (190) can thus be cast in the equivalent form 



C 



eff 



— n 



H — n k 
2 



(93) 



which turns out to be exact pZTQ. 

The theory describes a sound wave, with the dimensionless phase field <p rep- 
resenting the Goldstone mode of the spontaneously broken global U(l) symmetry. 
It has the gapless spectrum E 2 (k) = c 2 k 2 . The effective theory gives, ignoring 
vortices for the moment, a complete description of the superfluid at low energies 
and small momenta. When one goes to higher energies and momenta, additional 
terms with higher-order derivatives should be included in the effective theory, but 
it remains a phase-only theory. 



4.2 Coulomb potential 

It is straightforward to generalize the result (^UJ) to include long-ranged interac- 
tions. A case of particular interest to us is the 3-dimensional Coulomb potential 



V(x) = h, 



(94) 



whose Fourier transform in d space dimensions reads 



v(k) = 2 d -\ (d - 1)/2 riUd - 1)}- 



\d-l " 



(95) 



Here, q stands for the electric charge, which in the case of Cooper pairs is twice 
the electron charge. The simple contact interaction L mt = —\ B f x \ cj)(x)\ 4 in Eq. 
( f49l ) is now replaced by 



J 'mt 



1 

2 Jx,y 



|0(t,x)|V(x-y)|0(t,y) 



(96) 



The rationale for using the 3-dimensional Coulomb potential, even when consid- 
ering charges confined to move in a lower dimensional space, is that the electro- 
magnetic interaction remains 3-dimensional. The effective theory then becomes 
after passing over to the Fourier representation 



C cS = -nU(k) + ^U(k , k)y^jU(k , -k), 



and leads to the spectrum 

£ 2 (k) = 2 d -\ {d - 1)/2 T[l(d- 1)] 



n Q , k |3-d 

m 



(97) 



(98) 
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For d = 3, this yields the famous plasma mode, with an energy gap given by the 
plasma frequency cu 2 = Airnq 2 jm. For d = 2 on the other hand, the spectrum 

behaves as E(k) oc yjkj, implying that the mode it describes is much harder that 
the sound wave with the spectrum E(k) oc |k| obtained for the system without the 
1/r Coulomb interaction included. 

To appreciate under which circumstances the Coulomb interaction becomes 
important, we note that for electronic systems l/|x| ~ k F for dimensional rea- 
sons and the fermion number density n ~ fcf, where k-p is the Fermi momentum. 
The ratio of the Coulomb interaction energy to the Fermi energy eF = k\j2m is 
therefore proportional to n~ l l d . This means that the lower the electron number 
density is, the more important the Coulomb interaction becomes. 



4.3 Hyperscaling 

Let us consider the two terms in the effective theory ( ^3"] ) quadratic in the Gold- 
stone field ip and write them in the most general form [ ^8| ] 

•M--i4(V*,)» + i-» 
2 m z 2 



4? = -7~(Vv?) 2 + -n^{d^f. (99) 



The coefficient p s is the superfluid mass density, which is, as we saw in the pre- 
vious section, a response function and in general does not equal mn. The other 
coefficient, 

* 2 * = P L , (100) 

can be related to the (0,0)-component of the polarization tensor n o- This can 
be understood by noting that an electromagnetic field is included via the minimal 
substitution <9 M — > + qA^, with the electromagnetic vector potential. Since 
the polarization tensor (times q 2 ) is obtained by differentiating the effective theory 
twice with respect to the vector potential, we obtain 

lim n 00 (0,k) = n 2 n, (101) 

k| — >0 

where, as is typical for response functions, the energy transfer is put to zero before 
the momentum transfer k is. Equation (P9h leads to the general expression for the 



speed of sound 

m z n z K 

The singular behavior of the system close to the critical point is encoded in 
the phase-only theory. Simple dimensional analysis shows that near the phase 
transition it scales as 

Cos ~r (d+2) , (103) 
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while 

(vv?) 2 ~ r 2 , ov) 2 ~ c 2 ~ r 22 , d04) 

with £ t the correlation time and z the dynamic exponent. Combining these hyper- 
scaling arguments, and remembering that the mass parameter is inessential with 
regards to the critical behavior, one arrives [|28|] at the scaling laws for the two 
coefficients appearing in the effective theory ($2%: 

p s ~ ^-(«»+*-2) j K „ £-(«*-*). (105) 

The first conclusion is consistent with the universal jump predicted by Nelson and 
Kosterlitz [E9p which corresponds to taking z = and d = 2. 



In the presence of impurities it is believed that the compressibility stays finite 
at the critical point, implying z = d [Eft]. This remarkably simple argument thus 



predicts an exact and nontrivial value for the dynamic exponent. 

Without impurities, the dynamic exponent is z = 2 [|T7|]. This agrees with 
what one naively expects, given that in the nonrelativistic theory ( ft9| ) we started 
with, one time derivative appears in combination with two space derivatives, ido + 
V 2 /2m. This last argument should, however, be treated with care when applied to 
the phase-only theory (|9"3|). In that theory, the time and space derivatives appear in 
a symmetrical form, yet z is in general not unity, as we just saw. The difference is 
that in the effective theory, the relative coefficient c 2 scales according to Eq. ( |102| ) 
with the scaling laws ( jl05| ) as 

c 2 ~e 2(1 " z) , (106) 

while the relative coefficient m in the microscopic theory does not scale. Inciden- 
tally, the (quantum) XY model has a dynamic exponent z = 1, so that c in that 
case does not scale. 

In experiments on charged systems, instead of the superfluid mass density, 
usually the conductivity a is measured. To see the relation between the two, we 
introduce a vector potential in the effective theory by replacing Vy? with Vy2 — 
qA in Eq. (^9|), and allow the superfluid mass density to vary in space and time. 
The term in the action quadratic in A then becomes after passing to the Fourier 
representation 

S ° = I M-k)Ps(k)A(k). (107) 

2 m z Jk 

The electromagnetic current, 



j « - ^R) <108) 

obtained from this action can be written as 

= a(k)E(k), (109) 
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with the conductivity 



a(k) 



q 2 p s {k) 



(110) 



m? k 

essentially given by the superfluid mass density. So if we know the scaling of the 
electric charge, we can determine the scaling of the conductivity. 

With the 1/r Coulomb potential included, the quadratic terms in the effective 
theory (|97|) may, after passing to the Fourier representation, be cast in the general 
form 

1 / A, . o Ikl^ 1 



,(2) 
-cff 



Ps k 2 




2 q 
where q' is the redefined charge parameter 

q >2 = 2 d-l 7r (d-l)/2 r 

The charge is connected to the (0, 0)-component of the polarization tensor via 

i k id-i 



(HI) 



q 



(112) 



q 



a 



lim 



(113) 



IkHo n 00 (o,k) 

A simple hyperscaling argument like the one given above for the case without 



Coulomb interaction shows that near the transition, the charge scales as [ pOQ 

e~ z , (H4) 



J2 tl-z 

q ~ 



independent of the number of space dimensions d. It then follows from Eq. ( |1 1(J| ) 
that the conductivity scales as 



3~(d+z) 



(115) 



Exercise: Give an alternative derivation of the result ( 114), using Eq. (1101 



In the presence of random impurities, the charge is expected to be finite at 
the transition, so that z = 1 [j30|]. This is again an exact result, which replaces 
the value z = d of an impure system without Coulomb interaction. The pre- 
diction was first confirmed for impure superconducting films [Q], and has subse- 
quently also been observed in other 2-dimensional systems such 2-dimensional 
Josephson-junction arrays [J3T1], quantum Hall systems [|32"|], and 2-dimensional 



electron systems Q33| ] . We will refer to a quantum critical point with a 1/r Coulomb 
interaction as CQCP. In the vicinity of such a critical point, the conductivity scales 



as [34] 



2-d 



(116) 



implying that in two space dimensions, the conductivity is a marginal operator 
which remains finite at the CQCP. 
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4.4 Scaling of magnetic vector potential 

Let us finish this section by determining the scaling of the magnetic vector po- 
tential. We start with the observation that close to a CQCP, the electric field E 
scales as E ~ Ci 1 ^ 1 ~ (for a review, see Ref. Q350). Thus conductiv- 



ity measurements P5| , |3"3| ] close to a CQCP collapse onto a single curve when 
plotted as function of the dimensionless combination where as before 

S = K—K c ) measures the distance from the critical point K c , and v is the correla- 
tion length exponent, £ ~ \&\~ v ■ (For a field-controlled transition, K stands for the 
applied magnetic field, while for a density-controlled transition it stands for the 
charge-carrier density.) The scaling of the electric field with the correlation length 
expresses the more fundamental result that the anomalous scaling dimension d A 
of the magnetic vector potential A is unity, d A = 1. 

Because the magnetic vector potential always appears in the gauge-invariant 
combination V — qA, the anomalous scaling dimension of the electric charge q 
of the charge carriers times the vector potential is unity too, d qA = 1. Writing 
the anomalous scaling dimension of the vector potential as a sum d A = d A + 
|r^A of its canonical scaling dimension d% = \{d + z — 2), obtained by simple 
power counting, and (half) the critical exponent r/ A , describing how the correlation 
function decays at the critical point, we conclude that d q = d q — \r]x- Here, 
dg = 1 — d A stands for the canonical scaling dimension of the electric charge. 
Now, for a 1/r Coulomb potential, the charge scales according to Eq. (|1 14]) as 



q 2 ~ independent of the number d of space dimensions []30t] . Combined with 
the previous result, this fixes the value of the exponent r/ A in terms of the number 
of space dimensions and the dynamic exponent: 

VA = 5-d-2z. (117) 

Since in the presence of impurities, the electric charge is finite at a CQCP, leading 
to z = 1, it follows that i] A = 1 in two space dimensions. 

As we shall see in the next section, this exponent becomes important when 
considering the interaction between vortices close to the CQCP. 

4.5 Experimental status 

For a critical discussion of the experimental status of the phase-only theory, see 



Ref. Q37Q. A more recent discussion can be found in Ref. 

According to the phase-only theory discussed here, no electronic excitations 
exist in the critical region. However, electron tunneling measurements on super- 
conducting films of varying thickness apparently probed the energy gap of these 



excitations [39]. Moreover, the gap was found to approach zero as the transition 



to the insulating state is approached. Similar experiments [ |4"0| ] for the field-tuned 
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transition showed the presence of a large number of electronic excitations near the 
Fermi energy, thus raising doubts about the applicability of the phase-only theory. 

Experimental support for the presence of electron pairs in the insulating state 
comes from Hall effect studies on superconducting films, which show two critical 
fields [f|IJ]. The lower critical field is seen in the longitudinal resistance and is 
believed to mark the superconductor-insulator transition. The higher critical field 
is seen in the transverse or Hall resistance and is believed to signal the crossover 
from a bosonic to a fermionic insulator without pairing. At the higher critical 
field, the longitudinal resistance has its maximum. 

The critical exponents determined in earlier experiments on the superconductor- 
insulator transition [|2|, ^6)] had the value z = 1 for the dynamic exponent, in accord 
with the prediction in Ref. ATI], and v = 1.3 for the correlation length exponent. 
More recent studies p2|], however, find agreement with these results only for the 
transition tuned by changing the film thickness. For the field-tuned transition the 
value zv = 0.7 was found instead, which is about half the value one expects. The 
cause for this discrepancy is not clear. It implies that, contrary to common believe, 
the critical exponents depend on how the phase transition is crossed, by tuning the 
field or the film thickness. 

Clearly, more experimental and theoretical studies are required to fully under- 
stand the superconductor-insulator transition, and to establish to what extend the 
phase-only theory is applicable. 



5 Duality 

One of the most intriguing results found in experiments on quantum phase transi- 
tions in superconducting films, as well as in 2-dimensional Josephson-junction ar- 
rays [|3T|], quantum Hall systems [|43|], and 2-dimensional electron systems [ |33| ] is 



the striking similarity in the current- voltage (I-V) characteristics on both sides of 
the transition. By interchanging the I and V axes in one phase, an I-V characteris- 
tic of that phase at a given value of the applied magnetic field (in superconducting 
films, 2-dimensional Josephson-junction arrays, and quantum Hall systems) or 
charge carrier density (in 2-dimensional electron systems) can be mapped onto an 
I-V characteristic of the other phase at a different value of the magnetic field or 
charge-carrier density. This reflection symmetry hints at a deep connection be- 
tween the conduction mechanisms in the two phases that can be understood by 
invoking a duality transformation [[]], |4j]. Whereas the conducting phase is most 
succinctly described in terms of charge carriers of the system, the insulating phase 
is best formulated in terms of vortices. At zero temperature, these topological de- 
fects should, just like the charge carriers, be thought of as quantum point particles. 
The duality transformation links the two descriptions, which turn out to be very 



27 



similar. 



5.1 Vortices 

Let us now include vortices in the phase-only theory. This is achieved by introduc- 
ing the so-called plastic field cp p via the minimal substitution d^cp — > d^p + ip p 
p5[]. The plastic field is defined such that its curl gives a delta function at the 
location of the vortices. Specifically, in two space dimensions, where vortices are 
point objects, located at the positions x a say: 

V x cp p = -27r]T<5(x-x a ), (118) 

a 

while in three dimensions, where vortices are line objects, located along the curves 
C a say: 

V x ip p = -2nJ2 [ dx Q 5(x-x a ). (119) 

Let us concentrate on static phenomena so that we can ignore the time deriva- 
tives in the effective theory (|99|). When besides vortices also the magnetic vector 
potential is included, the effective theory becomes in three dimensions 



4? = - ^ - qA)2 - \ {v x A)2 ' (l20) 

or after the canonical transformation qA — ► qA — (p p : 

= -\%i<y<p - - \<y >< A - BP ) 2 ' ( 121 ) 

where the plastic field B p stands for 

B P = -$oE/ dx a 5(x-x Q ), (122) 

q, J Cry 

with $ = 2-k /q the magnetic flux quantum in units where the speed of light and 
Planck's constant % are set to unity. [In two dimensions, this plastic field is a 
scalar and stands for 

5 p = -<& £S(x-x a ), (123) 

a 

as follows from Eq. ( |1 18| ).] 

After integrating out the phase field tp in Eq. ( |121| ), we obtain the magnetic 
part of the effective action S mSLg . Written as a functional integral over the magnetic 
vector potential, it reads in the Coulomb gauge V ■ A = 



e i5mag = / exp 

A 



■-(Vx A-B p ) 2 -ii-A 2 

2 ^ 2 A 2 



(124) 
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with A the magnetic penetration depth. The mass term, with A~ 2 = q 2 p s /m 2 , is 
generated through the Anderson-Higgs mechanism in the process of integrating 
out the phase mode cp. 

With this construction, we can now calculate the interaction between two vor- 
tices. To facilitate the calculation in the case of a superconducting film below, 
we linearize the first term in Eq. ( |124| ) by introducing an auxiliary field h via a 
Hubbard-Stratonovich transformation: 

— |(V x A - B p ) 2 -> i(V x A - B p ) • h - |h 2 . (125) 

The original form is regained after integrating out the auxiliary field again. After 
integrating out the magnetic vector potential, we arrive at a form appropriate for a 
dual description in terms of magnetic vortices rather than electric charges [ ^6] ] 

e iS mag = J exp j f J [-IA 2 (V x h) 2 - |h 2 - th ■ B p ] | . (126) 

Physically, h, which satisfies the condition V • h = 0, represents (i times) the 
fluctuating local induction. The vortices described by B p couple to h with a 
coupling constant g = $o/A independent of the electric charge. Observe the close 
similarity between the original ( |124| ) and the dual form fll26p . This becomes even 
more so when an external electric current j p is coupled to the A field by including 
a term —A • j p in Eq. ( ]124[ ), and B p describing the vortices is set to zero there. 



Finally, also integrating out the local induction, one obtains the well-known 
Biot-Savart law for the interaction potential S mag = — J t V between two static 
vortices in a bulk superconductor [p7|], 

V(r) = ^£ yJ Bf(x)G(x-y)Sf(y) 
= 9 -l f dl^dl 2 ^ 



An Jci Jc 2 R 
= -|!L[ln(r/2A)+7]+0(r/A) 2 , (127) 

Z7T 

where we ignored the self-interaction. In Eq. ( |127j ), G(x) is the vortex-vortex 
correlation function with Fourier transform G(k) = l/(k 2 + A~ 2 ), R denotes 
the distance between the differential lengths dl 1 and dl 2 , L is the length of each 
of the two vortices, and 7 is Euler's constant. For distances smaller than the 
magnetic penetration depth, which is the length scale for variations in the electric 
current and the magnetic field, the interaction is logarithmic as in a superfluid. If 
the system size is smaller than A, it will replace the penetration depth as infra- 
red cutoff in the logarithm, and there will be no reference to the electric charge 
anymore. 
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To describe magnetic vortices in a film of thickness w [|4a|], the bulk result 
( |126[ ) has to be adjusted in two ways to account for the fact that both the vortices 
and the screening currents, which produce the second term in ( |126| ), are confined 
to the plane. This is achieved by including a Dirac delta function wS(x 3 ) in the 



second and third term. Instead of Eq. fll26Q , we then arrive at the interaction 
potential [|8|, [47]] 



V±(r) = [ B p ± (x ± )G ± (x ± -y ± )B p ± (y ± ) 

= _^i[ln(r/4A ± )+ 7 ] + 0(r/A ± ) 2 , (128) 

where B± = — $ J2 a $( x ± ~~ x ±) describes the vortices in the film with coordi- 
nates x_i_, \j_ = \ 2 /w is the transverse magnetic penetration depth, g\ = &l/\± 
the coupling constant squared, and 



G_l(x_l) = / G_ L (x._ L ,x 3 ) 

JX3 

= [ e~ ik ^G ± (k ± ,0), (129) 

Jk ± 

with (jrj_(kj_, 0) = 2/k±(2k± + Xj 1 ). For small distances, the interaction is seen 



to be identical to that in a bulk superconductor p8Q, and also to that in a superfluid 
film. As in the bulk, the vortex coupling constant g± in the film is independent of 
the electric charge. 

The logarithmic interaction between vortices we found in Eq. ( |128| ) appears to 
pose a severe problem to the duality picture we alluded to in the introduction of 
this section as the charges interact via a 1/r Coulomb potential. The difference 
should spoil the experimentally observed reflection symmetry in the I-V charac- 
teristics. However, it should be realized that the results derived in this subsection 
are valid only in the mean-field region, where T)a — 0. In the critical region gov- 



erned by a CQCP, the value of this exponent was found in Sec. to be unity. 
As we will now demonstrate, this leads to a qualitative change in the interaction 
potential between two vortices from logarithmic in the mean-field region to 1/r 
in the vicinity of the CQCP. 

5.2 Changing vortex interaction 

Close to a CQCP we have to include the field renormalization factor Za in the 
vortex-vortex correlation function G± appearing in the expression ( |128| ) for the 
vortex interaction. It then becomes 



2 Z 



G ± (k ± ,0) = ^ (130) 
k± 2k± + A_l 
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with Zx ~ k 1 ^ . Because the magnetic vector potential and the local induction 
renormalize in the same way, their renormalization factor is identical. Due to this 
extra factor, the interaction between two vortices in the film takes the form of a 
1/r Coulomb potential p^ ] 

VAr) = f -, (131) 
2tt r 

where a is some microscopic length scale which accompanies the renormalization 



factor Z A for dimensional reasons [p0|]. 

The absence of any reference to the electric charge in the renormalized and 
bare interaction (at least for small enough systems) implies that the same results 
should be derivable from our starting theory ( |120D with q set to zero. By directly 
integrating out the Anderson-Bogoliubov mode, and ignoring the momentum de- 
pendence of p s , which is valid outside the critical region, one easily reproduces 
the bare interaction potential fl!28p . The renormalized interaction ( |13ip is ob- 
tained by realizing that according to Eq. ( |105[ ), p s ~ k± for d = 2 and z — 1. In 
other words, the extra factor of k± that came in via the renormalization factor Z A 
in our first calculation to produce the 1/r potential, comes in via p s here. 

One might wonder if perhaps also the Coulomb interaction between electric 
charges changes in the vicinity of a CQCP. We do not expect this to happen. Since 
the 1/r Coulomb interaction is genuine 3-dimensional, this interaction cannot be 
affected too much by what happens in the film, which constitutes a mere slice of 3- 
dimensional space. The reason that the interaction between vortices is susceptible 
to the presence of a CQCP, is that this interaction is a result of currents around the 
vortex cores which are confined to the plane. 

A similar change in the r-dependence of the interaction between two vortices 
upon entering a critical region has been observed numerically in the 3-dimensional 
Ginzburg-Landau model [f5TJ] . Near the charged fixed point of that theory, r/ a = 1 



[p2Q , as in our case. 

This is a very pleasing coincidence as the (2+l)-dimensional Ginzburg-Landau 
model constitutes the dual formulation of the system. 

5.3 Dual Theory 

The dynamics of the charged degrees of freedom is described by the effective 
Lagrangian with the speed of sound given by Eq. ( |102[ ). In accord with the 



above findings, we have ignored the coupling to the magnetic vector potential, so 
that the effective theory essentially describes a superfluid. 

In the dual formulation, the roles of charges and vortices are interchanged. 
And the Anderson-Bogoliubov mode mediating the interaction between two vor- 
tices is represented as a photon associated with a fictitious gauge field a^, i.e., 
dpip ~ t^xdydx- In 2+1 dimensions, this identification makes sense as a photon 
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has only one transverse direction and thus only one degree of freedom. It there- 
fore represents the same number of degrees of freedom as does the Anderson- 
Bogoliubov mode. 

The elementary excitations of the dual theory are the vortices, described by a 
complex scalar field ip. Specifically, the dual theory of Eq. d99| ) turns out to be the 
Ginzburg-Landau model 



£duai = -\fl„ + \(d, - tgaM 2 - m%\^\ 2 - u\^\\ (132) 

with f^ u = d^a u — d u a^, a mass parameter, and u the strength of the self- 
coupling. Both the gauge part as well as the matter part of the dual theory are of 
a relativistic form. The gauge part is because the effective theory (p9|), obtained 
after ignoring nonlinear terms, is Lorentz invariant. The matter part is because 
vortices of positive and negative circulation can annihilate, and can also be cre- 
ated. In this sense they behave as relativistic particles. As was pointed out in Ref. 



[[44]], the speed of "light" in the gauge and matter part need not to be identical and 
will in general differ. 

The interaction potential ( |128| ) between two external vortices is now being 
interpreted as the 2-dimensional Coulomb potential between charges. The ob- 
servation concerning the critical behavior of the Ginzburg-Landau model implies 
that the qualitative change in V(r) upon entering the critical region is properly 
represented in the dual formulation. 

Whereas in the conducting phase, the charges are condensed, in the insulating 
phase, the vortices are condensed [O. In the dual theory, the vortex condensate is 
represented by a nonzero expectation value of the ip field, which in turn leads via 
the Anderson-Higgs mechanism to a mass term for the gauge field a^. Because 
( e fj,v\d u a\) 2 ~ (d^ip) 2 , the mass term a 2 ^ with two derivatives less implies that 
the Anderson-Bogoliubov mode has acquired an energy gap. That is to say, the 
phase where the vortices are condensed is incompressible and indeed an insulator. 
Since electric charges are seen by the dual theory as flux quanta, they are expelled 
from the system as long as the dual theory is in the Meissner state. Above the 
critical field h = Vj_ x a = h Cl they start penetrating the system and form an 
Abrikosov lattice. In the original formulation, this corresponds to a Wigner crystal 
of the charges. Finally, when more charges are added and the dual field reaches 
the critical value h C2 , the lattice melts and the charges condense in the superfluid 
phase described by the effective theory 
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